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Exchange corrections in a low temperature plasma 
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We have studied the exchange corrections to linear electrostatic wave propagation in a plasma using a quan¬ 
tum kinetic formalism. Specifically we have considered the zero temperature limit. In order to simplify the 
calculations we have focused on the long wavelength limit, i.e. wavelengths much longer than the de Broglie 
wavelength. For the case of ion-acoustic waves we have calculated the exchange correction both to the damping 
rate and the real part of the frequency. For Langmuir waves the frequency shift due to exchange effects is found. 
Our results are compared with the frequency shifts deduced from commonly used exchange potentials which 
are computed from density functional theory. 
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I. INTRODUCTION 


Recently much work has been devoted to quantum plas¬ 
mas, see e.g. the books and review articles 111-0]. The re¬ 
search is motivated by an interest in e.g. quantum wells |5|, 
spintronics |0, plasmonics 10. laser plasma interaction |8], 
astrophysical applications |0 or general theory development 
||I3]. The theoretical descrmtions range from hydrodynamical 
equations (e.g. Refs. |0 0 fTH] ) to quantum kinetic models 
(e.g. Refs. 10 ESlil [Tj]] and field theoretical approaches 
(e.g. Ref. lUill L Most models include the physical effects of 
particle dispersion and Fermi pressure, and in some cases the 
magnetic dipole force and magnetization due to the electron 
spin lift [iZ [iM[l^ . An important effect that is sometimes 
accounted for (e.g. Refs. (Tl], UtI - I^ ') but often overlooked 
is the exchange effects resulting from the total antisymme¬ 
try of the electron wavefunction. A popular approach to in¬ 
clude the effects of exchange interaction has been to apply 
density dependent potentials deduced from density functional 
theory (DFT) ifTHIl^l - l^ . An advantage with this is that the 
resulting fluid models becomes comparatively simple once the 
exchange potentials are established. As a consequence prob¬ 
lems involving both high-frequency dynamics j 11, [H as well 
as low-frequency (ion-acoustic) dynamics IMl ^21 can be ad¬ 


dressed in a straightforward way, also for nonlinear problems 
in. A drawback is that the calculation of DFT potentials 
typically involve approximations (e.g. the local density ap¬ 
proximation (FDA)) whose accuracy can be hard to estimate 
beforehand. Thus there is a general need to validate results 
derived from DFT by independent methods. 

In the present paper we calculate the exchange contribu¬ 
tion to the ion-acoustic dispersion relation using quantum ki¬ 
netic theory derived from first principles. Previous works 
along this line 110. [isl] have assumed the ordering T 3> 7 f 
( where T is the temperature and 7> is the Fermi temperature), 
which has prevented a direct comparison with results based 
on DFT potentials that have considered the opposite ordering. 
In this paper we focus on the low-temperature limit T ^ 7> 
and evaluate the exchange contribution to the ion-acoustic dis¬ 
persion relation in the Hartree-Fock approximation to first or¬ 
der in perturbation theory. We deduce that the effects of the 
exchange term is to increase the phase-velocity of the ion- 


acoustic mode and to increase the linear damping rate (which 
is due to wave-particle interaction). Moreover, as a confirma¬ 
tion of the correctness of the quantum kinetic formalism, we 
compare results from our quantum kinetic formalism with pre¬ 
vious results for the exchange contribution of high-frequency 
Langmuir waves. In this latter case we recover the results 
of Refs. 120 - 1^ exactly within the outlined approximation 
scheme. 

Finally we compare our findings with results based on com¬ 
monly used DFT potentials fillll9[ 1^ 1^ . As the DFT po¬ 
tentials are incorporated in a fluid formalism, no comparison 
can be made for the damping due to wave-particle interaction. 
In general we find a qualitative agreement. In particular the 
frequency shift in the different formalisms has the same scal¬ 
ing with the parameters (density and wavenumber) both for 
ion-acoustic and Langmuir waves. However, there is a dis¬ 
crepancy concerning a numerical factor. This is discussed in 
more detail in the final section of the manuscript. 


II, THE EXCHANGE CORRECTION AT T = 0 K 


In a previous paper. Ref. IItI] . the exchange contribution to 
the evolution equation of the Wigner function was derived (see 
Eq. (7) of Ref. HI). The correction was obtained by writing 
down the first equation in the BBGKY-hierarchy and writing 
the two-particle density matrix as a anti-symmetric product of 
one-particle density matrices. The treatment was here limited 
to electrostatic fields. For a generalization allowing for elec¬ 
tromagnetic fields, see Ref. S. Eq. (7) of Ref. d was 
further simplified by considering a plasma without spin polar¬ 
ization and summing over all spin states, and also by taking 
the long scale limit (where the macroscopic scale length is 
assumed to be much longer than the de Broglie wavelength). 
The long scale assumption implies that the Wigner function 
reduces to the Vlasov limit, in which case the Wigner func¬ 
tion becomes similar to a classical distribution function. This 
means that the evolution equation reduces to the Vlasov equa¬ 
tion with a correction term due to the exchange effects. The 
resulting expression (Eq. (12) of Ref. HtIi I with the exchange 
correction written in the right hand side reads 
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Here / denotes the electron distribution function, E is the 
electrical field, —e and nie is the electron charge and mass re¬ 
spectively, h = 2nh is Planck’s constant and we use x and r for 
position vectors and p and q for momentum vectors. Further¬ 
more, = d I dxi and analogously for dp and d^.. An arrow 
above an operator indicates in which direction it acts. We have 
also used the summation convention so that a sum over in¬ 


dices occurring twice in a term is understood. Finally V (r) = 
e^/47reo|r| is the Coulomb potential. Here we will use Eq. 
([T]i to consider linear wave propagation treating the exchange 
term on the right hand side perturbatively. Thus we linearize 
Eq. ([T]) and make a plane wave ansatz E^ = E exp(—/of + ikz), 
/(x,P,f) =/o(p)+/iexp(-/of-f/fe). 

To first order in a long wavelength expansion, the linearized 
version of Eq. O then reduces to 
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Previously the correction to the dispersion relations for ion 
acoustic waves (see Ref. 113 ]) and for Langmuir waves (see 
Ref. ifisll ) has been found considering a Maxwellian back¬ 
ground, by making the approximation fi « /i in the exchange 
term. The symbol f\ = qE /[/(o — kp^/me)\dfo/d denotes 
the solution when exchange effects are neglected. 

Here we will study the same problems but for the case when 
T = OK. As we will demonstrate the problem is analytically 
tractable owing to the simple form of the background distri¬ 
bution function /o which is now given by 

r.tp) _ / 2/(2;r/i)^ IPI<PF (3) 

0, |P|>PF, 


is the combined ion and electron susceptibility in the absence 
of exchange effects. Here we have neglected a contribution 
from a small but finite ion Eermi velocity (i.e. we have used 
that vpi <C (o/k). In the expression above COe and cOi denote the 
electron and ion plasma frequency respectively. The exchange 
correction is given by 


Zi = - 


with 
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where pF = is the Eermi momentum (we will 

also use the notation vf = Pf /me for the electron Eermi veloc¬ 
ity below), and iiq is the equilibrium electron number density. 
In particular, taking /i = /i in the exchange term, combining 
Eqs. (2) and @ with Poisson’s equation and computing the 
ion charge density classically we obtain 

1 =Zo + Zi, 

where 
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By changing integration variables to Ui = p + q, U 2 = p — q, 
with the Jacobian equal to 1/8, we can use the properties 
of the Dirac delta functions to reduce the problem to an in¬ 
tegral over two spheres. Introducing spherical coordinates 
u; = Ui (cos (pi sin 6i , sin <pi sin 0,, cos 9i ) we explicitly get Xi = 
m^q^I' /8;r®eok^/i^with I' given by 
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where y/ is the angle between ui and U 2 and a = (o/{kvp) cosines cosy/ = cos0icos02-f sin0i sin02cos(^i — ^). By 
is the dimensionless phase velocity. By the spherical law of 
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changing variables to <pi, ^ = <pi — <p 2 , the (pi integral is trivial. 
The (p integral can then be performed, giving a factor 7t j{a^ — 
where a = I — cos 0i cos 02 and b = sin 0i sin 02. Since 
a^ — b^ = (cos 01 — cos 02 )^, we have that/' = where 

I is the double integral 


/ = 
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At this point, no approximations have been made beyond 
those leading to the evolution equation (|2|i, taking /i = /i, 
and r = OK. For these approximations the exchange correc¬ 
tion is therefore proportional to (l9]l for all frequency regimes. 
Thus Poisson’s equation, with exchange corrections, is 
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where Xt) is given by Eq. (0). From now on we must treat 
the low-frequency ion-acoustic case separately from the high- 
frequency Langmuir case in order to simplify the expression 
for/. 

a. Ion-acoustic waves The integral / can be evaluated 
analytically in terms of a with computer algebra software. We 
will consider the quasi-neutral limit (co ^ 0),) where the left 
hand side of Eq. ( fTOl i is negligible. For ion-acoustic waves we 
can then make the approximation that a = when 

solving the integral (|9]l. It is also possible to plug in the value 
of a and evaluate the integral numerically. The dispersion 
relation for ion-acoustic waves with exchange corrections is 
then 


(0^ = 
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To the best of our knowledge this result has not been derived 
before. 

b. Langmuir waves In the high-frequency regime 
o/k 2> v'F, an expansion of I in powers of v'f k /(0 can be 
made. With the lowest order non-vanishing correction, the 
dispersion relation in this regime is 
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This is in exact agreement with previous calculations using 
several different methods, see Refs. ll24l - l^ . 


III. DFT COMPARISON 

Eq. (|9]l together with the specific results (fTTT i and (fT^ are 
the main results of the present paper. The exact agreement of 
(fT^ with previous results also confirms the validity of (|9]l. 
A strength of the quantum kinetic formalism is that it fol¬ 
lows from first principles, and that it can address wave par¬ 
ticle interaction, such as the enhanced Landau damping rate 
found in Eq. (dUl. However, a drawback is that the formal¬ 
ism is difficult to use for more complicated problems. Thus 


there is a need to develop theories that are easier to apply in 
a more general context. One such possibility is offered by 
DFT. A commonly used exchange potential (see e.g. Refs. 
llTll[l^l^l2^ ') computed from DFT is 

0.985(37r2)2/3;^2^2 /^^y/3 

4;r MeVp \noJ 

that has been used for low-frequency ion-acoustic phenomena 
(e.g. Refs, in, 1^1 ) as well as high-frequency Langmuir 
waves (e.g. Refs. Inllll]). Including the contribution from 
Eg. (foi l in the electron momentum equation (see e.g. Ref. 
11 nil I and treating the term as a small perturbation when calcu¬ 
lating the ion-acoustic dispersion relation, we get a qualitative 
agreement with Eq. (fTTl i. In particular the ion-acoustic fre¬ 
quency is decreased and this change scale as Tp'coj/m^Vp, in 
accordance with (fTTl i. The numerical value of the frequency 
shift deviates rather significantly from our result, however. We 
note that if we make the adjustment 0.985 —6.52 of the nu¬ 
merical pre-factor in ( fT3l l we would get agreement with the 
real part of the frequency in Eq. (fTTI) . Here we stress that the 
wave damping cannot be compared with the DFT formalism, 
as this requires a quantum kinetic framework. 

The same comparison can be made in the high-frequency 
regime. Thus we again include the contribution from ( fOl l in 
the electron momentum equation, calculating the frequency 
shift of the Langmuir dispersion relation. Also here we have 
qualitative agreement, i.e. the Langmuir frequency increases 
proportional to k^, and the scaling with density is in accor¬ 
dance with the factor lYcOg jm^Vp. This time the numerical 
accuracy is better, although not perfect, and we need to make a 
less significant substitution of the numerical factor in Eq. d, 
0.985 — ?> 1.23 in order to get agreement with our result (fT2l i. 
Thus we conclude that the DFT exchange potential Eq. (fT3l) 
are in qualitative agreement with the results obtained here, but 
that the numerical accuracy is better for Langmuir waves than 
for ion-acoustic waves. 


IV. SUMMARY AND CONCLUSION 

In the present paper we have computed the exchange contri¬ 
bution to the ion-acoustic dispersion relation in a plasma, us¬ 
ing a quantum kinetic formalism. The validity of our approach 
have been confirmed by comparison with similar results for 
high-frequency Langmuir waves 113,1251, in which case we 
get exact agreement for the exchange correction. This is to be 
expected, as our formalism as well as that of Refs, 
is based on first principles. While the quantum kinetic for¬ 
malism is of fundamental importance, it has the drawback of 
producing complicated formulas (cf. Eq. (1)) that can only be 
solved perturbatively. In fact even a perturbative treatment is 
far from straightforward. Thus there is a need for a formalism 
that can be used in more complex situations. Such a possibil¬ 
ity is offered by DFT, where the resulting exchange potentials 
can be used in fluid theories, and straightforward ly appli ed to 
a number of linear and nonlinear problems lUll ll^ 1^ 1^ . 
However, as the computation of DFT potentials typically in¬ 
volves approximations (e.g. the local density approximation 
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(LDA)) whose accuracy is unknown, there is a need to evalu¬ 
ate DFT potentials against independent methods. A key mo¬ 
tivation in the present pa per has been to evaluate the DFT 
potentials used in Refs, ill ll IT^ 1^ l23l against the quan¬ 
tum kinetic results computed in our formalism. This com¬ 
parison reveals that there is a reasonable qualitative agree¬ 
ment. In both cases the relative magnitude of the exchange 
term scales as /i^cUg/mgV^, and both the Langmuir and ion- 
acoustic wave frequencies decreases due to the exchange in¬ 
teraction, in agreement with Eqs. (fTTT i and (fT2l i. The numer¬ 
ical value of this frequency shift differ somewhat, however. 
To some extent this can be fixed by replacing the numerical 
pre-factor in the exchange potential. If this approach is cho¬ 
sen, we note that different substitutions must be used for low- 
frequency and high-frequency phenomena (i.e. 0.985 6.52 

in the former case, 0.985 —?> 1.23 in the latter case). It is an 
open question to what extent this result is robust (i.e. whether 
the same numerical coefficient is a good approximation for 
different problems at a given frequency scale), or if more ad¬ 
vanced expressions for the exchange potentials are needed to 


cover a broad spectrum of problems. 

The fact that there is not a perfect agreement between 
the DFT exchange potential and the quantum kinetic theory 
should not be overly surprising. The DFT potentials used in 
Refs. lUllll^l^B^ have not necessarily been optimized for 
the situation we have been studying. In our case the fields are 
dynamically varying (such that time-dependent density func¬ 
tional theory (TDDFT) applies), and the system is weakly col- 
lisional such that the collision-free Vlasov equation holds to 
leading order. The systems of relevance for the present study 
generally include plasmas with a high density and a modest 
temperature T <^Tf. In a laboratory context this applies to 
e.g. solid state plasmas and inertial confinement fusion plas¬ 
mas before the heating stage, and in an astrophysical context 
this include e.g. white dwarf stars. In that case the perturbed 
distribution function will be comparatively far from thermo¬ 
dynamic equilibrium, and approaches which do not include 
the full quantum kinetic details run the risk of losing informa¬ 
tion. Nevertheless an approach based on density functional 
theory is a valuable option, in particular for problems where 
the approach used here becomes too complicated. 
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